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ABSTRACT 

[LEC-671  linear  time  serial  algorithm  for  testing  planarity  of  graphs  uses  the  linear 
time  serial  algorithm  of  [ET-76]  for  .^-numbering.  This  ^-numbering  algorithm  is 
based  on  depth-first  search  (DFS).  A  known  conjecture  states  that  DFS,  which  is  a 
key  technique  in  designing  serial  algorithms,  is  not  amenable  to  poly-log  time 
parallelism  using  "around  linearly"  (or  even  polynomially)  many  processors.  The 
first  contribution  of  this  paper  is  a  general  method  for  searching  efficiently  in 
parallel  undirected  graphs,  called  ear-decomposition  search  (EDS). 
The  second  contribution  demonstrates  the  applicability  of  this  search  method.  We 
present  an  efficient  parallel  algorithm  for  5f-numbering  in  a  biconnected  graph.  The 
algorithm  runs  in  logarithmic  time  using  a  linear  number  of  processors  on  a 
concurrent-read  concurrent-write  iCRCW)  PRAM.  An  efficient  parallel  algorithm 
for  the  problem  did  not  exist  before.  The  problem  was  not  even  known  to  be  in 
NC. 

1.   Introduction 

We  define  the  problems  considered  in  this  paper.  For  all  these  problems  we  use  the 
same  input. 

Input.    An  undirected  graph  G(V,E)  and  some  specified  edge  e  =  (s,t)  in  E .  (Denote  «  =  |V| 
and  m  =  \E\.) 

Let  Pq  be  the  path  that  consists  of  the  edge  e .  An  ear  decomposition  of  G  starting  with 
PQ  is  a  decomposition  E  =  PQ  U  Pl  U  ■  •  U  Pk,  where  P,_t  is  a  simple  path  whose 
endpoints  belong  to  P0  U  •  ■  ■  U  P,,  but  its  internal  vertices  do  not.  A  simple  path  Pt  is 
called  an  ear.  It  is  called  an  open  ear  if  the  two  endpoints  of  P,  do  not  coincide,  and  a 
closed  ear  otherwise.  An  ear  decomposition  is  called  open  if  all  its  ears  are  open. 

The  ear  decomposition  problem.    Find  an  ear  decomposition  starting  with  P0. 

The  open  ear  decomposition  problem.    Find  an  open  ear  decomposition  starting  with  Pq. 

An  one-to-one  function  /  from  V  to  {1 /;}  is  called  an  st-  numbering  if  it  satisfies: 

(i)  f(s)=\    and  /(;)=«,   and  (ii)  for  each  viV— {s,t}  there  exist  adjacent  vertices  Vj  and  v2 
such  that/(v:X/(r)</(r;). 

The  st-numbering  problem.    Find  an  5f-numbering  of  G. 

A  graph  G  has  an  open  ear  decomposition  starting  with  an  edge  iff  G  is  biconnected 
(see  [W-32]).  A  graph  G  is  biconnected  iff  it  has  an  •rr-numbering  for  every  given  edge 
is .t)  ( see  [LEC-67]). 


The  model  of  parallel  computation  used  in  this  paper  is  the  concurrent-read 
concurrent  write  (CRCW)  parallel  random  access  machine  (PRAM).  A  PRAM  employs  p 
synchronous  processors  all  having  access  to  a  common  memory.  A  CRCW  PRAM  allows 
simultaneous  access  by  more  than  one  processor  to  the  same  memory  location  for  either 
reads  or  writes.  We  use  the  following  concurrent-write  convention.  In  case  several 
processors  seek  access  to  the  same  memory  location  for  write  purposes,  one  of  them 
succeeds  but  we  do  not  know  in  advance  which.  See  [V-83]  for  a  survey  of  results 
concerning  PRAMs. 

For  each  of  the  three  problems  mentioned  above  we  give  a  parallel  algorithm  which 
runs  in  0(\ogn)  time  using  n  +  m  processors  on  a  CRCW  PRAM.  Alternative  parallel 
implementations  attaining  optimal  speed  up  where  m>n\ogn  are  also  given  for  each 
problem.  These  alternative  algorithms  run  in  time  0(log«log("  «log( ~]n).  where  log1*'  is  the 
k'h  iterate  of  the  log  function. 

[L-85]  showed  that  the  ear  decomposition  problem  has  a  parallel  algorithm  which  runs 
in  poly-log  time  using -a  polynomial  number  of  processors  and  is  therefore  in  the  class  NC. 
A  remark  at  the  end  of  Section  2.1,  explains  why  we  believe  that  it  was  not  known  whether 
the  open  ear  decomposition  problem  is  in  the  class  NC.  We  are  also  not  aware  of  any 
parallel  algorithm  for  5f-numbenng  Apparently,  it  was  not  even  known  if  ?r-numbering  is 
in  NC.  Here,  we  do  not  only  determine  that  these  two  problem  are  in  NC,  but  actually  give 
very  efficient  parallel  algorithms  Each  of  these  algorithms  runs  in  logarithmic  time  using  a 
linear  number  of  processors. 

Serial  graph  algorithms  use  two  main  techniques  for  searching  graphs:  Depth-First 
Search  (DFS)  and  Breadth-First  Search  (BFSi.  Both  techniques  seem  to  be  not  ideal  for 
parallel  computation.  The  mo^t  efficient  poly-log  time  parallel  implementations  of  BFS 
require  a  number  of  processors  which  is  cubic  in  the  number  of  vertices.  Many  researchers 
have  conjectured  that  DFS  is  not  even  in  NC  One  of  the  most  challenging  tasks  in  parallel 
computation  is  to  cope  to  with  this  apparent  intractability  of  DFS.  One  approach  is  not  to 
give  up  and  try  to  get  the  most  out  of  DFS  as  in  [A-85]  This  approach  is  interesting  from 
the  theoretical  point  of  view  since  it  bounds  the  limits  of  DFS  in  parallel.  We  believe, 
however,  that  when  it  comes  to  actually  designing  efficient  parallel  algorithms,  a  more 
realistic  approach   should  be   taken     That   is.   BFS  and   DFS  should   be   replaced  b>   new 
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techniques  which  are  amenable  to  parallel  computation.  This  approach  of  looking  for  new 
techniques  to  replace  DFS  has  been  practiced  in  designing  efficient  parallel  algorithms  for 
biconnectivity  [TV-85]  and  strong  orientation  [V-85].  Observe,  however,  that  unlike  our 
problems,  it  was  trivial  to  establish  the  membership  of  biconnectivity  and  strong  orientation 
in  NC  using  simple  transitive  closure  algorithms. 

Ear-decomposition  has  the  flavor  of  a  general  search  technique  in  graphs.  It  arranges 
the  vertices  of  the  graph  by  partitioning  them  into  paths.  This  enables  further  exploration 
of  the  graph  in  an  "orderly"  manner  We  call  this  search  technique  ear  -  decomposition 
search  fEDS). 

For  more  motivation  for  considering  the  ear  decomposition  problems,  we  refer  the 
reader  to  [L-85]. 

The  serial  algorithm  for  <?-numbering  is  considered  as  one  of  the  classical  serial  graph 
algorithms  (see  [E-79]i.  The  linear  time  algorithm  for  the  problem  [ET-76J  is  heavily  based 
on  DFS  The  present  paper  copes  with  the  apparent  intractability  of  DFS  for  parallel 
computation  by  providing  a  new  algorithm  tor  the  problem.  The  new  algorithm  is  based 
on  the  new  EDS  technique  for  ear  decomposition  of  a  graph,  thereby  demonstrating  its 
applicability.  The  algorithm  is  also  based  on  refined  insights  into  the  .^-numbering 
problem.  Interestingly,  our  algorithm  provides  also  a  new  linear  time  serial  algorithm  for 
sr-numbering. 

The  sf-numbering  is  used  as  an  important  component  in  several  serial  algorithms. 
Most  known  is  the  planarity  testing  algorithm  of  [LEC-67]  which  was  improved  to  run  in 
linear  time,  following  the  linear  time  algorithm  for  .sr-numbering  of  [ET-76].  Our  algorithm 
gives  hope  for  finding  efficient  parallel  planarity  testing  using  the  approach  of  [LEC-67].  It 
is  yet  unresolved  whether  there  exists  an  efficient  parallel  algorithm  for  planarity  testing. 
Note  that  the  involved  planarity  testing  of  [JS-82]  still  leaves  much  to  be  desired,  as  it  runs 
on  O(\oz~n)  time  using  0(n~     /log";?  I  processors. 

Another  application  for  sr-numbering  is  mentioned  in  [IR-84]  for  the  following 
problem.  Given  an  undirected  biconnected  graph  fine  two  ^panning  trees  T  and  7;  which 
are  rooted  at  the  same  vertex  r  so  mat  the  paths  from  r  to  v  in  T  and  T~_  are  vertex 
disjoint  i  except  r  and  v)  This  problem  is  important  tor  reliable  communication  as  was 
indicated  in  that  paper 
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In  Section  2  we  present  a  parallel  algorithm  for  computing  an  ear-decomposition  and 
an  open  ear-decomposition.  In  Section  3  we  use  the  open  ear  decomposition  to  compute 
^/-numbering  of  a  biconnected  graph  G.  The  numbering  is  done  in  two  stages:  (1)  We 
orient  the  edges  of  the  input  graph.  As  a  result  we  get  a  directed  acyclic  graph  such  that 
each  of  its  vertices  is  on  a  path  from  5  to  r.  (2)  We  topo logically  sort  of  this  digraph  to  get 
the  .^-numbering.  The  whole  subtlety  of  the  algorithm  lies  in  the  orientation  of  Stage  1. 
Through  a  careful  case  analysis  we  perform  this  orientation  based  on  considerations  which 
are  essentially  local.  We  were  surprised  by  the  fact  that  these  local  considerations  were 
sufficient.  Another  intriguing  comment  is  that  one  of  the  interesting  open  problems  in 
parallel  computation  is  whether  there  is  a  poly-log  time  algorithm  for  topological  sort  that 
uses  linearly  many  processors.  We  have  examples  where  our  topological  sort  algorithm 
fails  on  some  acyclic  digraphs.  However,  we  can  prove  that  this  topological  sort  algorithm 
will  work  correctly  on  all  acyclic  digraphs  which  can  be  obtained  as  outcome  of  the  first 
stage. 

2.   Ear-Decomposition  In  Parallel 

Let  G(V,E)  be  an  undirected  biconnected  graph.  We  give  an  algorithm  for  finding  an 
open  ear  decomposition  of  G  starting  with  (s,t),  where  (s,t)dE.  A  parallel 
implementation  of  this  algorithm  is  described  at  the  end  of  the  section. 

Definitions:  Let  T(V.E)  be  a  tree  which  is  rooted  at  t. 

(1)  For  viV,  LEVEL(v)  is  the  length,  counting  edges,  of  the  unique  path  in  T  from  t  to  v. 

(2)  For  viV,  F(v)  is  the  father  of  v  in  T. 

(3)  For  u,v  i  V,  LCA(u,v)  is  the  lowest  common  ancestor  of  a  and  v  in  T. 

(4)  Let  u,v€V,  where    u  is  not  an  ancestor  of  v.  We  denote  by  euy  the  first  edge  in  the 
unique  path  in  T  from  LCA(u,v)  to  u. 

In  order  to  compute  an  open  ear-decomposition,  we  first  find  an  ear-decomposition  of 
G  which  is  not  necessarily  open.  This  is  described  in  the  first  subsection  below.  The  second 
subsection  shows  how  to  modify  this  ear-decomposition  into  an  open  ear-decomposition. 
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2.1.   Finding  an  ear-decomposition 

instead  of  designing  a  new  algorithm  for  finding  an  ear  decomposition,  we  use  a 
parallel  algorithm  for  a  different  problem  which  was  given  in  fV-85].  There,  the  algorithm 
assigns  directions  to  the  edges  of  a  connected  bndgeless  undirected  graph  so  that  the 
resulting  directed  graph  is  strongly  connected  We  note  that  we  obtain  the  same  ear- 
decomposition  as  in  [L-85].  However,  our  parallel  algorithm  is  more  efficient.  For 
completeness  of  the  presentation,  we  outline  below  the  algorithm  for  finding  an  ear- 
decomposition 

Step  t,  Find  a  spanning  tree  T(V.ET)  rooted  at  r  of  G .  such  that  (s.t)  will  be  the  only  tree 
edge  whose  endpoint  is  r.  This  is  done  as  follow^:  First,  find  a  spanning  tree  of  the 
subgraph  induced  by  the  vertices  V—  {t}.  Finally,  add  the  vertex  r  and  the  edge  (s.t)  to  the 
tree  'Remark:  The  fact  that  G  is  biconnected  implies  that  the  subgraph  induced  by  V-{r} 
is  connected.)  Step  1  is  the  only  step  in  this  subsection  which  is  not  completely  identical  to 
the  algorithm  of  [V-851.  There,  any  spanning  tree  will  do 

Edges  in  ET - {(.? ,r)}  will  be  referred  to  as  tree  edges  Edges  in  E  —  ET  will  be 
referred  to  as  non-tree  edges   The  edge  (s,t)  will  have  a  special  status. 

In  Step  2  we  assign  a  number  to  each  non-tree  edge. 

Step  2. 

(a)  For  each  vd  V.  compute  LEVEL(v)  and  F(v). 
For  each  non-tree  edge  («,v)  compute  LCAi  «,v). 

(b)  Assume  that  each  edge  etE  has  a  serial  number   1  <SERIAL(e)  <w      For  each  non- 
tree  edge  e ,  let  NUMBER {e)  =  ( LEVEL ( LCA ie\) SERIAL (e)). 

vV.  define  a  lexicographic  order  <L  on  these  numbers  as  follows  Let  e ,e'  be  non- 
tree  edges.  MUBERieXfSUUBER^e'  >  if  NUMBER(e).KNUMBER(e').l,  or 
NUMBER(e) .]  =  XUUBERi s'  •  !  and  SUMBERie  )  2<NUMBER(e').2,  where  NUMBER(x)A 
stands  tor  the  i'h  coordinate  of  the  pair  M'MBERi  x) 

..et  f  be  a  non-tree  edge  We  call  the  simple  cycle  which  is  formed  by  /  together  with 
the  edge*  in  Er  the  cycle  of  f 

Step  3.    Each  tree  edge  ,.     considers"  ail  non-tree  edge*  /  such  thai        is  in  their  cycle. 
Among  them  e  selects  an  edge  ;  whose  \'L'MBER(f)  is  minimal  (according  to  <L  ■■  to  be  its 
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master  edge.  The  edge  /  is  denoted  MASTER(e) . 

Proposition  2.1.1.  Each  non-tree  edge  e  together  with  all  the  tree  edges  which  selected  it  as 
their  master  edge  form  a  simple  path  or  a  simple  cycle.  We  call  this  path  or  cycle  the  ear 
of  e. 

Observe  that  the  order  <L  on  the  non-tree  edges  induces  an  order  on  the  defined  ears.  In 
addition,  define  the  edge  (s,t)  as  the  first  ear. 

Proposition  2.1.2.  The  defined  ears  and  the  order  yield  an  ear-decomposition. 

Note  that  some  of  these  ears  may  be  simple  cycles  and  therefore  the  resulted  ear- 
decomposition  may  be  not  open.  A  closed  ear  may  occur  when  a  non-tree  edge  (u,v)  is  a 
master  edge  of  all  the  tree  edges  in  its  cycle.  (For  an  example  where  an  arbitrary  order  on 
the  serial  numbers  of  the  non-tree  edges  imply  such  a  closed  ear,  see  edge  n2  in  Fig.  2.1. 
assuming  that  the  serial  number  of  n2  is  smaller  than  that  of  «1#) 

Remark:  Unlike  a  remark  in  Lovasz"s  paper.  Fig.  2.1  demonstrates  that  it  is  possible  to  get 
a  closed  ear  whose  endpoint  is  neither  the  root  nor  a  separating  vertex.  The  open  ear 
decomposition  proposed  in  [L-85]  is  based  on  this  erroneous  remark. 

2.2.    Finding  an  open  ear-decomposition 

In  order  to  have  a  decomposition  into  open  ears  we  have  to  define  a  subtler  order  on 
the  non-tree  edges.  We  show  how  to  do  it  by  reordering  the  non-tree  edges  which  have  the 
same  LCA.  Specifically,  this  reordering  is  achieved  by  changing  only  the  second  component 
in  the  pair  (LEVEL(LCAie))  .SERIAL!  e  )).  First,  we  take  a  closer  look  at  the  situation  where 
a  closed  ear  occurs.  Let  (u.v)  iE  —  ET  and  x  =  LCA(u.v)  and  suppose  xi=  u.  If  v  =  x  then  the 
ear  of  (u.v)  is  closed  iff  eux  chooses  (u.v)  as  its  master.  Otherwise  (i.e.,  v¥^x  and  hence 
evu  is  defined),  the  ear  of  («,v)  is  closed  iff  both  eu,  and  e,  u  choose  («,v)  as  their  master. 
Suppose,  for  instance,  that  euv  (or  e  u)  is  on  the  cycle  of  a  non-tree  edge  /  such  that 
LEVEL(LCA(f))<LEVEL(x).  Here,  no  matter  how  unfortunate  we  are  with  the 
initialization  of  SERIAL,  the  ear  of  e  cannot  be  closed  since  (u,v)  will  never  be  the  master 
of  eui.  Therefore,  we  can  characterize  all  the  difficult  cases,  as  the  ones  in  which 
LCA(MASTER(eu  ,,)  =  LCAi  MASTER!  e  ji  =  r  if  v*x  and  LCA(MASTER(eUA.))  =  x 
otherwise.  <MASTERt.\  is  the  ^ame  as  in  the  previous  subsection.)  This  leads  to  the 
following  observation,  which  motivated  us  in  changing  the  above  ear-decomposition  into  an 
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open  ear-decomposition. 

Observation.  Consider  any  assignment  of  numbers  into  SERIAL  of  the  non-tree  edges  so 
that  for  each  non-tree  edge  (m,v),  at  least  one  of  euv  or  ev  u  does  not  choose  (u,v)  as  its 
master.  Then,  the  above  ear-decomposition  must  be  open. 

This  change  is  done  using  an  auxiliary  (not  necessarily  connected)  graph  Hx  for  every 
vertex  r  in  V—  {t}. 

Step  4.  For  every  vertex  x  in  V—  {;},  we  construct  the  bipartite  undirected  graph 
Hx  =  (VX,EX)  ,  as  follows: 

Vx  =  {[u  ,v]  |(  u  ,v )  €  £  -  Er  and  LCA  ( u  ,v)  =  r}  (J  { [a  ,F(  u)]  \F( u )  =  x} . 

In  words,  the  vertices  Vx  are  all  non-tree  edges  whose  LCA  is  v  (to  be  called  non-tree 
vertices)  and  all  tree  edges  connecting  x  with  a  son  in  the  tree  (to  be  called  tree  vertices). 

Let  [«,v]  be  a  non-tree  vertex  in  Vx  such  that  u  is  not  an  ancestor  of  v  in  T.  Let  eux  be  the 
tree  edge  (w,je).  Then  the  edge  ([«,v],[w,jc])  is  in  Ex. 

The  graph  Hi  which  relates  to  the  graph  G  of  Fig.  2.1  is  given  in  Fig.  2.2. 

Step  5.  For  each  graph  Hx,  compute  its  connected  components  and  a  spanning  forest. 

Main  Lemma.  Consider  a  spanning  tree  of  a  connected  component  of  one  of  these  auxiliary 
graphs  Hx.  Then,  at  least  one  of  its  tree  vertices  [w,x]  (where  x  —  F(w))  must  satisfy 
LEVEL(LCA(MASTER(w,x)))<LEVEL(x). 

Relevance  of  Main  Lemma.  Let  [u,v]  be  a  non-tree  vertex  in  Hx  and  let  (w,x)  be  a  tree 
vertex  in  the  spanning  tree  of  the  connected  component  of  [u,v]  in  Hx  whose  existence  is 
guaranteed  by  the  Main  Lemma.  Consider  the  path  leading  from  [w,x]  to  [u,v]  in  this 
spanning    tree.  This  path  alternates  between  tree  and  non-tree  vertices.    Let  us  denote  it 

ti(  =  [w ,x]),nl,t2,n2 fa,«a(  =  [M,v]),  where  t,  stands  for  a  tree  vertex  and  n,  for  a  non-tree 

vertex.  Observe  that  the  edge  (represented  by)  ;,  belongs  to  the  cycle  of  (the  edge 
represented  by)  nh  for  l</<a,  and  to  the  cycle  of  n,_i,  for  2sf<a.  We  aim  that  none  of 
the  ears  of  n\,n2....,na  will  be  closed.  In  view  of  the  above  Observation,  we  simply  take 
care  that  the  edge  of  each  tree  vertex  t.  will  not  select  the  edge  of  the  non-tree  vertex  n;  as 
its  master.  For  this.  Step  6  below  numbers  the  non-tree  vertices  on  this  path  in  ascending 
order  (i.e..  so  that  SERIAL(n-J  <  <  SERIAHn  x)) .    This  implies  that  for  2<;^a.  the 
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of  t,  will  not  select  the  edge  of  «,  as  its  master.  In  addition,  the  Main  Lemma  implies  that 
LEVEL(LCA(MASTER(w,x)))<LEVEL(x),  implying  that  tx  which  belongs  to  the  cycle  of  n} 
is  not  included  in  the  ear  of  ri\. 

Proof  of  Main  Lemma.  Take  any  spanning  tree  as  in  the  Lemma.  The  proof  proceeds  as 
follows.  We  define  a  non  empty  subgraph  Gx  of  G  which  is  a  function  of  this  spanning 
tree.  We  show  that  if  the  Lemma  does  not  hold  with  respect  to  this  spanning  tree  then 
Gi^G  and  G]  is  a  biconnected  component  of  G.  This  would  contradict  the  biconnectivity 
of  G  and  therefore  the  Lemma  is  correct.  Let  us  define  G  X{V  y,E  \)  ■  Recall  that  all  the  tree 
vertices  in  the  spanning  tree  are  of  the  form  (v,F(v)),  where  F(v)=x.  The  vertices  of  Gj 
are:  (1)  All  the  vertices  v  such  that  [v,x]  is  a  tree  vertex  in  the  spanning  tree.  (2)  All 
vertices  of  G  which  are  descendants  of  the  vertices  in  Item  (1)  in  the  spanning  tree  of  G 
which  was  found  in  Step  1.    (3)  x  itself. 

G[  is  induced  by  these  vertices.  That  is,  the  edges  of  G\  are  all  the  edges  of  G  which  are 
incident  to  two  vertices  in  Gj. 

Suppose  Hx  has  more  than  one  connected  component.  Let  G^Vi^i)  be  a  similar  subgraph 
of  G  which  is  defined  for  the  spanning  tree  of  another  connected  component  of  Hx.  We 
first  show  that  it  is  impossible  to  have  an  edge  (u,v)  in  E  connecting  any  vertex  //  of 
V^-f.v}  with  any  vertex  v  of  V':-{.v}.  Clearly,  LCA(u,v)  must  be  x  itself.  The  edge  eu  v  is 
in  the  connected  component  of  Gy  in  Hx  and  evll  is  in  the  connected  component  of  G2  in 
Hx.  Therefore,  if  (u,v)  exists,  then  it  must  be  connected  to  the  tree  vertices  of  eu  ,  and  evu 
and  therefore  Gj  and  G?  are  the  same,  contrary  to  our  assumption. 

Hence,  if  there  is  no  edge  (u.v)  in  E  such  that  u  is  a  vertex  of  V,-{x}  and 
LEVEL(v)<LEVEL(x)  in  the  tree  of  G.  then  the  removal  of  x  will  disconnect  G.  This 
contradicts  the  biconnectivity  of  G  . 

Step  6. 

i  a)  For  the  spanning  tree  of  each  connected  component  in  each  auxiliary  graph  H..  find  a 
tree  edge  (w,x)  'where  F(w)=x)  such  that  LEVELiLCAi  MASTERl  u\.v) ) )  <LEVELi  v). 
(Its  existence  is  guaranteed  by  the  Main  Lemma.)  This  tree  edge  will  be  referred  to  as 
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the  edge  of  its  component 

(b)  Root  each  such  spanning  tree  at  the  edge  of  its  component  Compute  the  numbering  of 
the  non-tree  edges  in  a  preorder  traversal  of  the  spanning  tree.  For  every  non-tree 
edge  e,  this  numbering  will  be  stored  in  PREORDER(e) 

(c)  For  each  non-tree  edge  e,  let  NEWNUMBER(e)  =  (LEVEL(LCA(e)),PREORDER(e)). 

For  example,  notice  that  the  edge  of  the  component  given  in  Fig.  2.2  is  fj. 

Step  7.  Each  tree  edge  reselects  a  master  edge  according  to  XEWXUMBER  in  the  same  way 
as  the  selection  of  a  master  edge  in  Step  3  above. 

Proposition  2.2.1.  Each  non-tree  edge  e  together  with  all  the  tree  edges  which  selected  it  as 
their  master  edge  form  a  simple  path  (which  is  not  a  cycle).   This  path  is  the  ear  of  e. 

Again,  observe  that  the  order  <L  on  the  non-tree  edges  induces  an  order  on  these  ears.  In 
addition,  define  the  edge  (s.t)  as  the  first  ear. 

Proposition  2.2.2.  These  ears  and  order  yield  an  open  ear-decomposition 

Remarks: 

il)  This  algorithm  can  be  extended  to  check  if  G  is  biconnected.  In  case  that  G  is  not 
biconnected  it  can  be  modified  to  compute  the  open  ear-decomposition  of  each  biconnected 
component. 

(2)  For  the  ^-numbering  algorithm,  we  may  omit  each  non-tree  edge  e  of  G  which  forms 
an  ear  that  consists  of  e  only. 

2.3,    Parallel  Implementation 

We  first  present  an  implementation  which  runs  in  OdogMi  time  using  n+m 
processors.  Appendix  1  presents  an  alternative  implementation  which  runs  in  time 
dogf/log'  ^log'-'rt)  using  an  optimal  number  of  processors  when  m  >nlosn. 

Steps  1.  2  and  3  are  essential!;*  taken  from  [V-85].  These  steps  take  O(logn)  time 
using  n  +  m  processors. 

Remark.  The  algorithm  of  [V-85]  implements  steps  2  and  3  in  three  stages:  i  !  I  Fach  vertex 
v  selects  a  single  non-tree  edge  (u.v)  sucn  that  LCA(u.\  •  is  the  lowest  common  ancestor  of 

(v.v< .-..).  where  v;.  .        .vk  are  all  the  vertices  which  are  adjacent  to  v  by  an  edge  in 

E  —  Ej.    ■  2)  The  original  serial  numbers  of  the  non-tree  edges,  which  were  set  arbitrarily, 
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are  modified  in  the  following  way.  All  the  non-tree  edges  which  were  not  selected  in  Stage 
1  are  assigned  serial  numbers  which  are  greater  than  those  assigned  to  the  selected  non- 
tree  edges.  (3)  Each  tree  edge  selects  its  master  as  in  Step  3  above.  In  [V-85],  it  is  shown 
that  this,  rather  tedious,  implementation  is  more  efficient.  It  requires  only  n  computations 
of  LCA  instead  of  m  computations  as  in  the  direct  implementation. 

Finding  the  edges  of  the  Hx's  in  Step  4  can  be  done  in  O(losn)  time  using  n  +  m  processors. 
Note  that  the  total  number  of  edges  in  all  auxiliary  graphs  is  0(n  +  m).  Step  5  takes 
O(logn)  time  using  n  +  m  processors  using  the  parallel  connectivity  algorithm  of  [SV-82]. 
Step  6(a)  (i.e.,  finding  the  edge  of  each  component)  needs  0(1)  time  using  n  processors, 
since  we  check  for  each  tree  edge  whether  it  satisfies  the  inequality  of  Step  6(a).  Step  6(b) 
requires  rooting  of  each  spanning  tree  and  computation  of  preorder  numbering.  Using  the 
Euler  tour  technique  on  trees,  which  has  already  been  used  extensively  in  our  above 
references  to  [V-85],  this  can  be  done  in  0(\ogn)  time  using  n  +  m  processors. 

(  Remark:  Step  6(b)  requires  an  adjacency  list  representation  of  the  spanning  forest.  To 
obtain  adjacency  lists  from  a  list  of  /  edges,  one  may  consider  applications  of  the  sorting 
algorithm  of  [AKS-83]  that  uses  /  processors  to  sort  /  elements  in  time  <9(log/)  (when 
edges  are  compared  according  to  a  lexicographic  order).  Since  we  only  need  to  sort 
numbers  which  are  in  a  restricted  domain  we  can  apply  the  recent  randomized  radix  sort 
algorithm  of  [R-85]  which  achieves  almost-surely  the  same  asymptotic  time  (with  much 
more  acceptable  constants)  using  only  //log/  processors.  A  third  possibility  is  to  perform 
this  sorting  in  time  O(logZ)  and  /  processors,  using  an  adaptation  of  the  simple  notion  of 
"orthogonal  trees"  (see  [Th-83]).  However,  this  takes  more  space.    ) 

Step  6(c)  needs  0(\)  time  using  m  processors.  Step  7  can  be  implemented  in  the  same  way 
as  Step  3  above.  (Recall  that  for  the  ir-numbering  we  may  omit  the  edges  which  form  an 
ear  of  length  one.  This  can  be  done  in  0(log«)  time  using  m/\osn  processors  by  applying  a 
variation  of  a  parallel  "prefix  sums"  algorithm.  See  [FL-80].) 

3.   ST-Numbering  In  Parallel 

In    this    section    we    describe    how   to    compute    in    parallel    the    u-numbermg   of   a 
biconnected  graph  G.  given  an  open  ear-decomposition  of  G 
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The  input  of  our  algorithm  is  as  follows: 

1  li      A  graph  G(V,E)  and  a  specified  edge  (s,t)€E, 

(2)     An  open   ear  decomposition   G  =  Pq  U  P:  U  UP,  such   that  P0  =  (s,t)   and  the 

endpoints  of  each  P,  are  in  Pj  and  Pk,  where  ;>/>&>0.  (Observe  that  j  and  A-  depend 
^n  i>  j=j(0,k  =  k(i) .')  The  decomposition  is  given  in  the  following  form.  Each  edge 
in  Pj  knows  the  index  i,  and  has  pointers  to  its  neighbor  in  P,  leading  to  P  and  to  its 
neighbor  in  P,  leading  to  Pk.  Unless  otherwise  stated,  whenever  we  say  an  ear.  we 
will  mean  one  of  these  open  ears. 

Definitions: 

(  1)  The  endpoint  of  an  ear  P,.  which  belongs  to  P.  (resp.  Pk),  is  denoted  by  LiP,)  (resp. 
R(P,)).  (Recall  that,  in  the  above  notation,  j^k.)  In  addition,  s  (resp.  i)  is  defined 
to  be  L(Pq)  (resp.  P(P0))- 

(2)  The  vertex  in  an  ear  P,,  (si,  which  is  adjacent  to  HP.)  (resp.  P(P,)),  is  denoted  by 
LS(Pi)  (resp.  R5(P,-)). 

(3)  The  vertex  L(Pt)  is  called  the  anchor  of  P  . 

(4)  A  vertex  v  of  an  ear  P,  (denoted  by  v€P,)  will  be  called  an  internal  vertex  of  P,,  if  v 
is  neither  L(P,)  nor  R(P  ). 

Remark.  In  our  algorithm  for  ?r-numbering,  each  ear  is  "responsible  for  providing"  the 
numbering  of  its  internal  vertices  onh  Therefore,  we  can  ignore  ears  P  which  include 
only  one  edge  The  ear  Pn.  whose  only  edge  is  (s,l),  is  the  only  ear  of  this  type  which  is 
not  ignored.   This  means  that  we  have  O(n)  relevant  edges  in  the  graph  G . 

In  the  next  subsection  we  give  a  high-level  description  of  the  algorithm  A  description 
of  the  parallel  implementation  is  given  in  the  sequel. 

3.1.   High-level  description  of  the  algorithm 

The  numbering  is  done  in  two  stages.  In  Stage  1  we  orient  the  edges  of  the  input  graph 
so  that  the  resulting  directed  graph  is  acyclic  and  each  of  its  vertices  is  on  a  path  from  s  to 
r.  In  Stage  2  we  compute  a  topological  sort  of  this  digraph  This  results  in  assigning 
numbers  from  i  to  n  to  the  vertices  of  the  digraph  These  numners  yield  a  valid  st- 
numbering. 
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Remark:  The  topological  sort  of  Stage  2  is  applicable  to  digraphs  obtained  in  Stage   1,  but 
not  to  general  acyclic  digraphs. 

For  the  orientation  of  the  edges  of  Stage  1  we  use  the  ear-tree  ET(VT,ET)  which  is 
defined  as  follows: 

VT  =  {P!  |  P,  is  an  ear  of  G  }.    That  is,  every  ear  is  a  vertex  in  ET. 
ET  =  {(Pl,PJ)  \L(P,)  is  an  internal  vertex  of  Ps  ,  where  fSl}. 

Observation:  ET  is  a  directed  tree,  rooted  at  P0. 

The  descendance  relation  in  the  tree  ET  extends  naturally  to  the  vertices  of  G.  In 
particular,  if  an  ear  P,  is  a  descendant  of  an  ear  P;  (it  might  be  P,  itself)  in  ET.  we  say  that 
each  vertex  in  P,  is  a  descendant  of  L(Pi),  the  anchor  of  Pt.  Conversely,  such  L{Pt)  is 
called  an  anchor  — ancestor  of  each  vertex  in  P,. 

In  order  to  motivate  the  orientation  of  the  ears  in  Stage  1,  which  is  quite  involved,  we 
start  by  presenting  the  main  ideas  of  Stage  2  (the  numbering  algorithm).  This  will  be 
followed  by  a  description  of  Stage  1. 

The  principles  of  the  numbering  algorithm  (Stage  2). 

We  assume  that  each  ear  P,  O'^l)  was  oriented,  in  Stage  1,  either  from  HP,)  to  R(P,) 
or  from  R(P,)  to  LiP,).    The  ear  Pq  was  oriented  from  5  to  r. 

Definition:  An  ear  P,  which  is  directed  from  L(P,)  to  R(Pt)  (resp.  from  P(P,)  to  L{P,))  is 
called  an  outgoing  (resp.  incoming)  ear  of  L(P,). 

We  give  two  rules  which  imply  a  unique  numbering  of  the  vertices  of  G.  These  rules 
imply  actually  a  serial  algorithm  which  provides  this  numbering.  We  refer  often  in  our 
presentation  to  it  as  the  serial  numbering  algorithm.  Parallel  implementation  of  both  stages 
is  discussed  in  the  next  subsection. 

Rule  1  relates  to  the  relative  numbering  of  the  internal  vertices  of  an  ear. 

Definition:  Let  P,  be  an  ear.  and  let  u.  v  be  internal  vertices  of  P,.  We  say  that  u  precedes 
v  if.  when  moving  on  P,  according  to  its  direction,  we  encounter  u  before  encountering  v. 

Rule  1.  Let  P,  be  an  ear.  and  let  u,  v  be  internal  vertices  such  that  u  precedes  v.  The  st- 
numbering  of  u  and  each  of  its  descendants  is  smaller  than  the  sr-numbenng  of  v  and  each 
of  its  descendants. 
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Rule  I  guarantees  that  for  each  ear  P,,  /si,  all  its  interna!  vertices,  excluding  LS(P,) 
and  RS(P,i,  are  assigned  a  valid  ,?r-numbering.  This  holds,  since  each  of  these  vertices  will 
have  a  neighbor  in  P,  whose  number  is  smaller  and  a  neighbor  in  P,  whose  number  is 
greater  than  its  own  number.  Hence,  the  difficult)  in  providing  sr-numbering  is  reduced  to 
assigning  a  valid  vr-numbering  to  LS(P,)  and  RS(Pt),  for  each  ear  P,  Note  however,  that  if 
Pi  is  an  incoming  (resp.  outgoing)  ear  with  more  than  one  internal  vertex,  then  Rule  1 
guarantees  that  LStP^)  will  have  a  neighbor  whose  number  is  smaller  iresp.  greater)  than 
its  own  number,  and  RS(P,)  will  have  a  neighbor  whose  number  is  greater  (resp.  smaller) 
than  its  own  number. 

Rule  2  relates  to  the  relative    numbering  of  ears  which  have  the  same  anchor. 

Definition:  Let  v  be  an  internal  vertex  of  some  ear.  We  define  a  total  order  on  the 
incoming  and  outgoing  ears  of  v  as  follows: 

(i)  An  incoming  ear  /  of  v  is  before  another  incoming  ear  /'  of  v  iff  the  edge  (v.LS(l)) 

is  before  (v,LS(I'))  in  the  adjacency  list  of  v  in  G     (Any  other  arbitrary  order  on 
the  edge:,  will  do.  We  only  need  some  consistent  order.) 

(ii)         An   outgoing   ear   O  of  v    is   before   another   outgoing   ear   O'   of  v   iff  the   edge 
(v,LS(0))  appears  after  (\\LS(0'))  in  the  adjacency  list  of  i-  in  G 

(iii)        Each  incoming  ear  /  of  v  is  before  each  outgoing   ear  O  01  v. 

Remark:  Note  that  the  internal  order  defined  on  the  outgoing  ears  is  opposite  to  the  internal 
order  defined  on  the  incoming  ears.  Interestingly,  this  subtle  definition  is  essential  for  the 
correctness  of  the  algorithm. 

Rule  2.    Let  v  be  an  internal  vertex  of  some  ear 

Rule  2.1.    The  vr-numbering  assigned  to  each  of  the  vertices  belonging  to  incoming 
ears  of  v  is  smaller  than  the  ?r-numbenng  assigned  to  v. 

Ruk-  2.2,  The  ^-numbering  assigned  to  v  is  smaller  than  the  .^-numbering  assigned  to 
each  of  I   e  vertices  belonging  to  outgoing  ears  of  v. 

Rule  2.3.  If  an  ear  P  of  v  is  before  another  car  P'  of  then  the  fr-numbering  assigned  to 
each  ol  the  internal  vertices  in  P  is  smaller  than  the  rr-numbering  assigned  to  each  of  the 
internal  vertices  in  P' 
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Consider  now  a  serial  implementation  of  the  numbering  which  obeys  Rule  2.  Upon 
arriving  at  a  vertex  v,  we  first  assign  numbers  to  each  of  its  incoming  ears  (and  their 
descendants)  in  order.  Then  we  number  v,  and  finally  we  assign  numbers  to  each  of  its 
outgoing  ears  (and  their  descendants)  in  order.  (We  use  the  order   defined  above.) 

Rules  2.1  and  2.2  guarantee  that  for  each  incoming  (resp.  outgoing)  ear  P,  which 
have  more  than  one  internal  vertex.  LS(Pt)  is  assigned  a  valid  .u-numbering;  Since  the 
number  of  L(Pt)  is  greater  (resp.  smaller)  than  the  number  of  L5(P,);  Recall  that  Rule  1 
implies  that  the  number  of  the  neighbor  of  LS{P,)  in  P;  is  smaller  (resp.  greater)  than  the 
number  of  LS(P,).  Therefore,  the  difficulty  in  providing  sr-numbering  is  further  reduced 
to  assigning  valid  sr-numbering  to  RS(P,),  for  each  ear  P,.  This  will  be  the  task  of  the 
orientation  algorithm  (which  will  make  use  of  Rule  2.3). 

We  already  mentioned  that  Rules  1  and  2  define  a  unique  numbering  of  the  vertices  of 
the  oriented  graph.  This  serial  numbering  algorithm  starts  at  ear  P0,  which  is  directed 
from  s  to  r,  and  arrives  to  each  ear  through  the  ear-tree  ET.  An  ear  is  always  "entered" 
through  its  anchor. 

Fig.  3.1(a)  gives  an  example  of  a  graph  with  directions  on  the  ears.  Specification  of 
/?(.)  for  all  ears  is  omitted.  Fig.  3.1(b)  demonstrates  the  result  of  an  application  of  the 
numbering  process. 

The  orientation  of  the  ears  (Stage  1). 

We  should  remember  that  the  orientation  of  the  ears  (Stage  1),  which  is  described 
below,  is  aimed  at  the  numbering  algorithm  (Stage  2),  which  was  outlined  above.  As  was 
mentioned  before,  the  main  task  of  the  orientation  algorithm  is  to  guarantee  that,  for  each 
ear  P,,  the  number  which  will  be  assigned  to  P(P,)  is  consistent  with  the  number  which  will 
be  assigned  to  RS(Pt).  Specifically,  this  means  that:  (i)  If  P,  gets  a  direction  from  HP,)  to 
R(Pt)  then  the  task  is  to  take  care  that  the  number  of  P(P,)  will  be  greater  than  the  number 
of  RS(Pj).  (ii)  If  P,  gets  a  direction  from  P(P,)  to  LiP,)  then  the  task  is  to  take  care  that  the 
number  of  R(P,)  will  be  smaller  than  the  number  of  RSiP,). 

The  relation  between  the  numbers  assigned  to  P(P,)  and  RS(P.)  is  determined  by  the 
following  event:  Will  the  serial  numbering  algorithm  visit  the  (internal  vertices  of  the)  ear 
P,  before  or  after  visiting  the  vertex  PiP,)1  Our  crucial  observation  is  that  we  can  make 
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this  event  depend  only  on  the  direction  which  will  be  given  to  some  single  ear  Pm.  As  will 
become  clear  later,  Pm  is  always  closer  than  P,  to  P0  in  the  ear-tree.  tPm  is  not  necessarily 
an  ancestor  of  P,.)  This  ear  Pm  will  be  called  the  hinge  of  Pt.  (Observe  that  m  is  a  function 
of  i,  m  =  m(i).) 

Definition:  We  say  that  the  ear  Pt  gets  the  same  direction  as  its  hinge  Pm,  if  our  orientation 
algorithm  implies  that  both  are  directed  either  from  L( .)  to  /?(.)  or  from  R(.)  to  L(.). 
Similarly,  the  ear  P,  gets  a  direction  opposite  to  Pm,  if  the  algorithm  implies  that  one  of 
them  is  directed  from  L( .)  to  /?(.)  iff  the  other  is  directed  from  /?(.)  to  L( .) . 

Below  we  describe  how  the  hinge  Pm  is  found,  and  how  to  determine  whether  Pt  will 
have  the  same  or  the  opposite  direction  as  Pm.  Since  this  description  is  not  short  we  show 
first  Step  2  (the  last  step)  of  Stage  1.  That  is,  how  to  determine  the  orientation  of  each  ear, 
once  we  know  for  each  ear  its  hinge  and  its  direction  with  respect  to  the  hinge. 

We  define  the  hinge -tree  HT(VHT,EHT)  as  follows: 

Vht~  {Pi  I  Pi  is  an  ear  ot  G  }.   That  is.  every  ear  is  a  vertex  in  HT. 

EHT={(Pi,Pm)  i  Pm  is  the  hinge  of  P,} 

Fact:  HT  is  a  directed  tree  rooted  at  Pq. 

Step  2.  Let  P0  be  directed  from  L(P0)  to  R(P0)  (i.e.,  from  s  to  f).  For  each  /s  1.  consider 
the  ear  P,  and  the  path  from  P,  to  P0  in  HT.  Count  the  number  of  ears  on  this  path  which 
get  a  direction  opposite  to  their  hinge.  If  this  number  is  even  then  P,  is  directed  from  L(P,) 
to  R(P,),  otherwise,  it  is  directed  from  P(P,)  to  L(Pt). 

Step  1.  Finding  the  hinge  Pm  of  ear  P, 

Definition:  We  say  that  a  vertex  v€P,  appears  to  the  left  of  «€P,  on  P,  if  v  appears  prior  to 
u  when  moving  on  P,  from  Z.(P,)     to  R(P,). 

For  finding  the  hinge  of  P, ,  we  consider  several  cases. 

Case  1.  The  vertices  P(P,)  and  L(Pt)  are  both  on  the  same  ear  P.,  see  Fig.  3.2. 

Instructions:  In  this  case,  the  ear  P.  is  set  to  be  the  hinge  of  P,.  and  P  gets  the  same 
direction  as  P.  iff  L(Pt)  appears  to  the  left  of  P(P.)  in  P,. 
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Discussion:  Notice  that  we  use  the   fact  that  the  given  ear-decomposition  of  G   is  open, 
therefore,  R(Pt)*L(P,). 

Since  this  is  the  first  case  we  consider,  we  elaborate  on  it  slightly.  In  later  cases,  the 
justification  of  the  choice  of  the  hinge  and  of  the  relation  between  the  orientations  of  an  ear 
and  its  hinge  will  be  shorter. 

In  Fig.  3.2(a),  P,  gets  the  same  direction  as  P r  On  the  other  hand,  in  Fig.  3.2(b)  it 
gets  a  direction  opposite  to  P y  To  justify  this  orientation,  recall  the  numbering  algorithm. 
Assume  that  the  situation  is  as  in  Fig.  3.2(a).  The  ear  Pi  will  be  directed  from  L(P,)  to 
P(P,)  iff  Pj  will  be  directed  from  L(Py)  to  R(Pj).  Now  if  both  conditions  hold,  the  number 
assigned  to  RS(Pt)  will  be  smaller  than  that  of  R(P,),  and  if  none  holds  then  the  number 
assigned  to  RSiP,)  will  be  greater  than  that  of  R(P,).  This  yields  a  consistent  numbering  for 
RS(Pj).  Similarly,  assume  that  the  situation  is  as  in  Fig.  3.2(b).  The  ear  Pt  will  be  directed 
from  Z.(P,)  to  RIP,)  if f  Pj  will  be  directed  from  R(Pj)  to  HPj).  If  both  conditions  hold,  the 
number  assigned  to  RS(Pj)  will  be  smaller  than  that  of  P(P,),  and  if  none  holds  then  the 
number  assigned  to  RS(Pj)  will  be  greater  than  that  of  R(P,).  This  will  provide,  again,  a 
consistent  numbering  for  RS(P,). 

Case  2.  The  vertices  Z.(P,)  and  R(P,)  belong  to  different  ears. 

Assume  that  L(P,)  is  in  Pj  and  R(P,)  is  in  Pk,  and  let  the  lowest  common  ancestor  of 
Pj  and  Pk  in  the  ear-tree  be  Pa. 

Case  2.1.    Pa  is  neither  P }  nor  Pk,  see  Fig.  3.3. 

Let  \-j  (resp.  vk)  be  the  anchor-ancestor  of  P }  (resp.  Pk)  in  Pa. 

Case  2.1.1.    v,#vt  ,  see  Fig.  3.3(a). 

Instructions:  (Similar  to  Case  1.)  The  ear  Pa  is  set  to  be  the  hinge  of  P,,  and  P,  gets  the 
same  direction  as  P a  iff  v<  appears  to  the  left  of  vk  in  Pa. 

Establishing  the  validity  of  the  numbering  of  RS(P,)  in  this  case  is  similar  to  Case  1. 

Case  2.1.2.    \>j  =  vk.  see  Fig.  3.3(b). 

Instructions:  Consider  the  first  edge  (Prt.P:  )  on  the  path  in  ET  from  Pa  to  P..  This  edge 
corresponds  to  an  edge  of  G  such  that  one  of  its  endpoints  is  v  .  Denote  this  edge  by  e,.  In 
the  same  way.  Pk-  and  ek  are  defined.  Since  v.  =  vk:=v,  both  edges  appear  in  the  adjacency 
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list  of  v  in  G.  If  e;  appears  before  ek  in  this  list,  then  the  hinge  of  P,  is  set  to  be  P.-  and  P, 
gets  a  direction  opposite  to  Py.  If  ek  appears  before  er  then  the  hinge  of  P,  is  set  to  be  Pk> 
and  Pj  gets  the  same  direction  as  Pk. 

Discussion:  The  following  claim  shows  that  also  in  this  case  the  vertex  RSiPt)  is  assigned  a 
valid  ir-numbering. 

Clam  I .  The  ear  P,  is  directed  from  L(P.)  to  P(P,)  (i.e.,  P,  is  an  outgoing  ear)  iff  the 
number  assigned  to  RS(P,)  is  smaller  than  the  number  assigned  to  P(P,) 

To  prove  Claim  1  we  need  the  total  order  defined  on  the  incoming  and  outgoing  ears 
of  a  vertex  v.  This  is  since  this  order  determines  the  order  of  numbering  the  vertices  in  the 
serial  numbering  algorithm  (by  Rule  2.3.)  The  crucial  point  is  that  Claim  1  holds 
regardless  of  the  directions  which  will  actually  be  given  to  P,   and  Pk-. 

Proof  of  Claim  I.    We  have  two  possibilities  to  consider. 

Possibilin-  (a):  The  edge  e.  appears  before  the  edge  t\.  in  the  adjacency  list  of  v. 

In  Possibility  (a),  P,  and  each  of  its  descendants  are  numbered  before  Pk  and  each  of  its 
descendants  iff  P,  is  an  incoming  ear  (regardless  of  the  direction  of  Pk).  But.  the 
instructions  of  Case  2.1.2  imply  that  P  is  an  incoming  ear  iff  P,  is  an  outgoing  ear. 
Recall  now  that  RS(P,\  is  a  descendant  of  P  .  and  R(Pj)  is  a  descendant  of  Pk-.  Combining 
the  above  together  we  obtain  that  P,  is  an  outgoing  ear  iff  the  number  assigned  to  RSiP,) 
is  smaller  than  the  number  assigned  to  R(P,),  as  required. 

Possibility  (b):  The  edge  ek  appears  before  the  edge  e,  in  the  adjacency  list  of  v. 

In  Possibility  (b),  P}  and  each  of  its  descendants  are  numbered  before  Pk-  and  each  of  its 
descendants  iff  /%•  is  an  outgoing  ear  (regardless  of  the  direction  of  Py).  But.  the 
instructions  of  Case  2.1.2.  imply  that  Pk-  is  an  outgoing  ear  iff  Pl  is  an  outgoing  ear. 
Recall,  again,  that  RS(Pt)  is  a  descendant  of  P.  .  and  R(P.)  is  a  descendant  of  Pk  . 
Combining  the  above  together  we  obtain  that  P,  is  an  outgoing  ear  iff  the  number  assigned 
to  RS(  Pj)    is  smaller  than  the  number  assigned  to  R(  P, ) ,  completing  the  proof  of  the  claim. 

Case  2.2.  P^  is  P..  .  ?ee  rig.  .;  -  (N'ote  that  P,x  can  not  be  P.  as  this  implies  that  j<k. 
while  the  definition  of  L\.P  i  and  RP>  implies  that  j^lc.) 
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Let  \  be  the  anchor-ancestor  of  P,  in  Pa.  Using  the  same  arguments  as  above,  we  get 
the  following. 

Case  2.2.1.    v#/?(P,),  see  Fig.  3.4(a). 

Instructions:  In  this  case  the  hinge  of  P,  is  set  to  be  Pa,  and  P,  gets  the  same  direction  as  Pa 
iff  v  appears  to  the  left  of  R(P,)  in  Pa. 

Establishing  the  validity  of  the  numbering  of  RS(P,)  in  this  case  is  similar  to  Case  1. 

Case  2.2.2.    \=R(Pl),  see  Fig.  3.4(b). 

Instructions:  Let  P y  be  the  first  ear  on  the  path  in  ET  from  Pa  to  Pj.  The  hinge  of  P,  is  set 
to  be  P    ,  and  P,  gets  a  direction  opposite  to  P }  . 

Discussion:  The  ear  P,  will  be  directed  from  L(Pt)  to  R(Pt)  iff  P}-  will  be  directed  from 
R(Pj-)  to  L(Pj').  If  both  conditions  hold,  the  number  assigned  to  RS(P,)  will  be  smaller 
than  that  of  R(P,){  =  Z.(P,<))>  and  if  none  holds  then  the  number  assigned  to  RS(Pt)  will  be 
greater  than  that  of  R(P,)(  =  HP ,■)) .  This  will  provide  a  consistent  numbering  for  RS(P,). 

This  completes  Step  1  of  Stage  1. 

The  above  presentation  shows  that  all  the  vertices,  excluding  s  and  f.  will  be  given  a 
valid  .sr-numbering.  In  order  to  establish  correctness  of  the  algorithm  it  remains  to  show 
that  s  and  t  are  also  assigned  a  valid  .sr-numbering. 

Claim  2.  The  algorithm  assigns  the  number  1  to  the  vertex  5,  and  the  number  n  to  the 
vertex  t. 

Proof  of  Claim  2.  Recall  that  P0  's  directed  from  s  to  t.  Consider  an  ear  P,  which  is  a  son 
of  P0  in  the  ear-tree.  The  definition  of  L(Pt)  and  R(P,)  implies  that  both  L(P,)  and  R(Pt) 
must  belong  to  P0.  (This  definition  states  that  if  L(P,)  is  in  P}  and  R(P,)  is  in  Pk  then 
;>A\)  This  means  that  all  the  ears  Pt  such  that  L(Pt)  =  s  (resp.  L(P,)  =  t)  must  satisfy 
RiP.)  =  i  (resp.  R(P ,)  =  <;).  Thus,  by  Case  1  of  the  orientation  algorithm,  for  all  P,  such  that 
LiP,)  =  s.  Pt  is  an  outgoing  ear,  and  for  all  P,  such  that  L(P,)  =  t,  Pt  is  an  incoming  ear. 
Hence,  from  Rules  2.1  and  2.2.  it  follows  that  1  is  assigned  to  5.  and  n  is  assigned  to  t. 
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3.2.    Parallel  Implementation 

We  describe  a  parallel  implementation  of  the  two  stages  of  the  .^-numbering 
algorithm. 

Stage  1.  In  Stage  1  we  start  by  constructing  the  ear  tree  ET  Let  /  be  the  number  of  ears 
that  have  at  least  one  internal  vertex.  Observe  the  /<«.  Finding  the  edges  of  ET  from  the 
output  of  the  open  ear  decomposition  algorithm  can  be  done  in  Oil)  time  using  n 
processors.  Computing  the  adjacency  lists  of  ET  needs  0(log/n  time  with  n  processors 
using  integer  sorting  algorithms,  e.g.  [AKS-83;.  References  to  other  possible  sorting 
algorithms  may  be  found  in  the  remark  at  the  end  of  Section  2.  Finding  the  hinges  requires 
/  computations  of  lowest  common  ancestor  in  ET.  This  can  be  done  in  OMog/i  time  using  / 
processors,  see  [V-85].  The  other  computations  of  Stage  1.  including  the  computation  of  the 
direction  of  the  ears  based  on  the  hinge  tree,  can  also  be  done  in  0(log/)  time  using  / 
processors.   Therefore,  Stage  1  takes  0(\ogn)  time  using  n  processors. 

Stage  2.  We  now  describe  how  to  implement  the  numbering  rules  in  parallel.  We 
construct  a  tree,  whose  preorder  numbering  obeys  the  numbering  rules.  This 
numbering -tree  XT  is  defined  as  follows.  Let  v  be  a  vertex  of  G.  The  vertex  v  induces 
two  vertices  vin  and  vou,  in  NT.  Also.  vin  and  vout  are  both  connected  by  an  edge.  The  rest 
of  the  edges  of  AT  are  defined  by  the  following  rules: 

(1)  Assume  that  v  is  in  P,  and  let  u  be  the  neighbor  of  v  in  P,  which  precedes  v  (if  such  a 
vertex  exists).  Then  uout  is  connected  to  vin. 

(2)  For  each  incoming  ear  Pt  whose  anchor  is  v  the  vertex  r.„  is  connected  to  n-,-n,  where 
w  =  RS(P,). 

For  each  outgoing  ear  P,  whose  anchor  is  v  the  vertex  vout  is  connected  to  u-m,  where 
w  =  LS(Pi). 

Fig.  3.5  illustrates  the  construction  of  AT  for  an  example  of  a  part  of  an  ear-tree. 

We  now  perform  a  preorder  traversal  of  \T  starting  at  the  root  sin.  We  observe  that 
the  order  m  which  vertices  of  the  form  vou,  are  visited  during  the  traversal  obeys  the  rules 
of  the  numbering  aigcruhm.  That  is.  tor  each  v  in  G  we  first  visit  v!n  and  all  the  vertices  ol 
incoming  ears  whose  anchor  is  \  and  their  descendants.  Then  we  visit  .,.,,.  Finally,  we  visit 
ail  the  vertices  of  outgoing  ears  whose  anchor  is  v  and  their  descendants.    The  order  in 
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which  the  incoming  (and  outgoing)  ears  of  v  are  visited    is  as  defined  in  Rule  2.    Note  that 
for  each  vertex  v  in  G,  its  vertices  in  AT  are  traversed  in  the  order:  first  v,„,  and  then  vout. 

0;ring  the  traversal  we  assign  preorder  numbers  to  the  vertices  of  the  form  voul, 
ignoring  all  the  vertices  of  the  form  v,„.  These  numbers  provide  the  sr-numbenng  of  the 
graph  G 

The  .onstruction  of  AT  can  be  done  in  0(\ogn)  time  using  n  processors.  The  preorder 
numbering  can  be  done  using  the  Euler  tour  technique  for  computing  tree  functions  given 
in  [V-85]  in  Oy\ogn)  time  using  n  processors. 

To  sura  up,  given  a  biconnected  graph  G.  we  present  a  parallel  algorithm  which  finds 
an  jr-numbering  of  G.  The  first  part  of  the  algorithm  finds  an  open  ear-decomposition  of 
G  This  is  done  either  in  O(\ogn>  time  using  m+n  processors,  or  in  OOognlog'-'/ilog'-'ni 
rime  using  ((nlog«  +  m)/\ogn\og{l)n\og{^)n)  processors.  The  second  part  finds  an  st- 
numbenng  of  G  in  Oi\ogn\  time  using  n  processors.  Therefore,  the  whole  .rr-numbering 
algorithm  takes  eicher  0(\ogn)  time  using  n  +  m  processors  or  0(log^?logl2)«log(31«)  time 
using  ((nlogn  +/«  )/\6gn\ogl'-)n[og[~'ln)  processors. 
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APPENDIX  1.  An  Alternative  Parallel  Implementation  of  the  Open 
Ear-Decomposition  Algorithm. 

Observe  that  the  total  number  of  operations  in  the  parallel  implementation  of  Section 
2  is  Oinlogn  +  mlogn),  while  the  number  of  operations  in  the  respective  serial  algorithm  is 
0(n  +  m).  In  the  following,  we  present  an  Oi\ogn\og{'^n\og{  'n)  time  implementation  in 
which  the  total  number  of  operations  is  only  0(n\ogn  +  m).  This  is  done  using  the  new 
algorithms  for  connectivity  and  linked-list  ranking  of  [CV-86].  Notice  that  Steps  1,4,5  and 
6  of  Section  2  are  the  only  steps  where  0(m\ogn)  operations  are  performed.  (Recall  the 
Remark  concerning  Steps  2  and  3  in  Subsection  2.3.) 

Step  I.  Step  1  can  be  implemented  in  Odognlog1  'rtlog'  '«)  time  with 
((n\ogn  +  m  )/lognlogl''/ilog(  'n)  processors  using  the  connectivity  algorithm  of  [CV-86]. 

Step  4.  Recall  that  in  Step  4  we  construct  the  auxiliary  graphs.  Note  that  for  doing  this  we 
have  to  compute  LCA(u,v),  euv  and  e,  u  for  each  non-tree  edge  (u,v).  Specifically,  for  each 
non-tree  edge  (u,v)  it  takes  O(\ogn)  time  using  one  processor  to  compute  LCA(u,v).  This 
implies  a  total  of  0(m\ognt  operations.  We  want  to  implement  this  step  in  only 
0(n\ogn  +  m\  operations.  We  modify  the  auxiliary  graphs  of  Subsection  2.2.  This  enables 
dispensing  with  explicit  LCA  computations  for  most  of  the  non-tree  edges  in  constructing 
the  auxiliary  graphs. 

For  every  vertex  x  in  V—{t},  we  construct  the  bipartite  undirected  graph  Hx  =  (VX,EX)  in 
the  following  way. 

The  vertices  in  Vx  are  of  three  types: 

(1)  All  tree  edges  connecting  x  with  a  son  in  the  tree  (the  tree  vertices). 

(2)  (Let  u  be  a  vertex  in  G.  We  define  Ad-LCA(u)  to  be  the  lowest  common  ancestor  in  T 

of  («,«!,  .  .  .  ,U[),  where  u^ u,  are  all  the  vertices  which  are  adjacent  to  u  by  an  edge 

in  E~Ej.)  All  non-tree  edges  (u,v)  such  that  LCA(  u.v)  —x  and  x  is  also  either  Ad-LCA(u) 
or  Ad-LCA(v).    (These  are  the  non-tree  vertices.) 

(  3)  A  new  vertex  x.  I  We  call  it  the  hyper—  vertex.) 

Let  [u.v]  be  a  non-tree  vertex  in  Vx  ^uch  that  u  is  not  an  ancestor  of  v  in  T.  It  v  is 
Ad—LCACtt)  then  the  edge  ([u.v},[w,x])  !  where  [w,x]  is  euj  is  in  £,.  Otherwise,  [i.e.. 
when  x  *  Ad-  LCAi  a)  .  the  edge  l[u.v],x)  is  in  Ex.  (Note  that  in  this  case  Ad-LCA(y)  is  a 
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descendant  of  Ad-LCA(u).) 

For  these  modified  auxiliary  graphs  we  have  the  following  Lemma 

Modified  Main  Lemma.    Consider  a  spanning  tree  of  a  connected  component  of  one  of  the 

auxiliary  graphs  Hr.   Then  either   it  contains  x  or  at  least  one  ol  its  tree  vertices  [w,x] 

(where  x  =  F(w))  must  satisfy  LEVEL{LCA{MASTER(w,x)))<LEVEL{x). 

Notice  that  if  the  edge  ([u,v],x)  is  in  Ex  then  the  ear  of  (u,v)  must  he  open. 

We  describe  how  to  construct  these  auxiliary  graphs  in  parallel.  First  we  compute 
Ad— LCA(v)  for  each  vertex  v  This  is  done  as  in  the  implementation  of  Step  2  (given  in 
[V-85])  in  Odogm  using  n  processors.  For  each  vertex  v.  we  also  compute 
Sub-Ad-LCA(v)  which  is  the  first  vertex  on  the  path  from  Ad-LCA(v)  to  v.  This  also  can 
be  done  in  O(\ogn)  time  using  n  processors. 

For  a  given  non-tree  edge  (u,v),  if  [u.v]  is  a  vertex  of  an  auxiliary  graph  Hx  then  x  has 
to  be  either  Ad-LCA(u)  or  Ad-LCA(v).  This  gives  at  most  two  auxiliary  graphs  in  which 
[«,v]  can  be  a  vertex.  It  takes  0(1)  time  using  a  single  processor  to  identify  these  two 
auxiliary  graphs. 

Suppose  x  =  Ad-LCA(u).  The  vertex  [u,v]  is  in  Ht  only  if  x  =  LCA(u,v).  We  find  out 
whether  x  =  LCA(u,v)  by  checking  whether  v  is  a  descendant  of  Sub  -Ad- LCA(u).  [V-85] 
shows  that  given  two  vertices  u  and  i  in  G.  it  is  possible  to  determine  whether  v  is  a 
descendant  of  u  in  F  in  0(1)  time  using  one  processor.  (This  is  done  using  the  preorder 
and  postorder  numbering  of  the  vertices,  which  were  computed  in  the  implementation  of 
Step  2  in  [V-85].)  Thus,  we  can  determine  whether  [u.v]  is  in  Hx  in  0(1)  time  using  one 
processor.  Suppose  x  =  Ad~LCA{u)  and  [u,v]  is  a  vertex  in  Hx.  We  show  how  to  find  the 
the  adjacent  vertices  of  [w,v]  in  Hx.  Recall  that,  by  the  definition  of  the  auxiliary  graphs. 
[«,v]  is  connected  to  the  tree  vertex  [Sub-Ad-LCA(u),x].  Suppose  further  that  v  is  not  an 
ancestor  of  u  (i.e.,  x±v).  If  <  =  Ad-LCA(v)  then  [«,v]  is  connected  by  an  edge  to  the  tree 
vertex  [Sub-Ad-LCA(v),x],  otherwise  [u,v]  is  connected  to  x.  So.  finding  the  adjacent 
vertices  of  [u,v]  in  HK  takes  constant  time  using  one  processor. 

The  case  x  =  Ad—LCA(v]  is  treated  similarly 

The  parallel  connectivity  algorithm  of  [CV-861  needs  an  adjacency  list  representation 
ot  the  auxiliary  graphs  Hx.  We  show  how  to  obtain  such  representation.  Note  that  naive 
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sorting  of  the  edges  can  not  be  applied  within  the  complexity  bounds  that  we  seek. 

The  adjacency  lists  of  the  non-tree  vertices  can  be  obtained,  as  described  above,  in  constant 
time  using  one  processor. 

We  show  how  to  compute  the  adjacency  lists  of  the  tree  vertices.  Let  v  be  a  vertex  in  G 
and  let  x  be  Ad-LCA(v).  Recall  that  a  non-tree  edge  (u,v)  is  a  non-tree  vertex  in  Hx  if 
and  only  if  x  =  LCA(u,v).  Moreover,  if  [u,v]  is  a  non-tree  vertex  in  Hx  then  it  is  connected 
by  an  edge  to  the  tree  vertex  [Sub— Ad— LCA(v),x].  So,  for  each  vertex  v  in  G  we  consider 
all  edges  (k,v)  such  that  LCA(u,v)  =  Ad—LCA(v)  as  one  block  and  apply  a  sorting 
algorithm  to  these  blocks  for  obtaining  the  adjacency  lists  of  the  tree  vertices.  This  sorting 
takes  0(\osn)  time  using  n  processors  since  we  have  only  n  blocks.  Thinning  out  these 
blocks  to  include  only  relevant  edges  is  easy  and  can  be  done  by  a  parallel  prefix  sum 
computation  (given  the  adjacency  list  representation  of  G)  in  O(logn)  time  using  m/logn 
processors. 

The  computation  of  the  adjacency  list  of  the  hyper-vertices  is  similar.  Let  v  be  a  vertex  in 
G  and  let  x  be  Ad-LCA(v).  Recall  that  a  non-tree  edge  (u.v)  is  connected  by  an  edge  to  x 
in  Hx  if  and  only  if  x  =  LCA(u.v)  ^Ad-LCA(u).  So.  for  each  vertex  v  in  G  we  consider  all 
edges  (u.v)  such  that  LCA(u,v)  =  Ad-LCA(v)¥=Ad-LCA(u)  as  one  block  and  proceed  as 
for  the  tree  vertices. 

Step  5.  In  Step  5  we  compute  the  connected  components  and  a  spanning  forest  of  each 
auxiliary  graph.  Note  that  the  total  number  of  vertices  in  these  graphs  is  0{m).  Thus,  it 
seems  that  O(m)  operations  will  not  be  enough.  However,  we  can  perform  only  O(m) 
operations  using  the  following  Fact. 

Fact.    The  degree  of  each  non-tree  vertex  in  the  auxiliary  graphs  is  at  most  two. 

Thus,  we  may  shrink  each  non-tree  vertex  [u,v]  of  degree  two  and  the  two  edges 
emanating  from  it  into  a  new  edge  connecting  the  two  neighbors  of  [u.v].  The  resulting 
graphs  have  O(n)  vertices  and  O(m)  edges,  and  hence  we  may  apply  the  connectivity 
algorithm  of  [CV-86]  to  get  the  same  results  as  in  Step  1.  Notice  that  the  connected 
components  and  the  spanning  forest  of  the  original  graphs  can  be  recovered  from  the  the 
connected  components  and  the  spanning  forest  of  the  shrunken  graphs  in  constant  time 
using  Olm)  processors. 
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Step  6.  This  step  can  be  implemented  using  the  Euler  tour  technique  applying  the  new 
optimal  ranking  algorithm  of  [CV-86]  in  time  6>(log«log(:)«)  using  m/\ogn\og{2)n 
processors. 

In  summary,  we  have  shown  an  implementation  which  runs  in  0(\ogn\og(2)nlogl2)n) 
time  with  ((nlogn  +  m)/\ogn  log(2,«log'3)/i)  processors. 
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